Spontaneous demagnetization of a dipolar spinor Bose gas at ultra-low magnetic field 
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Quantum degenerate Bose gases with an internal degree of freedom, known as spinor condensates, 
are natural candidates to study the interplay between magnetism and superfluidity. In the spinor 
condensates made of alkali atoms studied so far, the spinor properties are set by contact interac- 
tions, while magnetization is dynamically frozen, due to small magnetic dipole-dipole interactions. 
Here, we study the spinor properties of S=3 52 Cr atoms, in which relatively strong dipole-dipole 
interactions allow changes in magnetization. We observe a phase transition between a ferromagnetic 
phase and an unpolarized phase when the magnetic field is quenched to an extremely low value, 
below which interactions overwhelm the linear Zeeman effect. The BEC magnetization changes due 
to magnetic dipole-dipole interactions that set the dynamics. Our work opens up the experimental 
study of quantum magnetism with free magnetization using ultra-cold atoms. 

PACS numbers: 03.75.-b, 67.85.-d, 03.75.Mn 
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Within optical dipole traps, it is possible to trap all 
Zeeman states of an atomic species in its electronic 
ground state. This has enabled the field of multi- 
component (spinor) Bose- Einstein condensates JBECs), 
i.e. BECs with internal degrees of freedom [HQ. With 
short range interactions, collisions between atoms are 
described by various scattering lengths, which leads to 
spin-dependent contact interactions. New phases arise, 
which have been investigated for F = 1 and F = 2 
atoms, using Rb and Na, by studying miscibility 0], spin 
dynamics [3-6], and spin textures [7|. In these experi- 
ments, spin-dependent contact interactions are extremely 
small compared to the linear Zeeman effect. In addition, 
the gas magnetization remains constant for all practi- 
cal purposes because contact interactions are isotropic, 
and anisotropic dipole-dipole interactions between alkali 
atoms are negligible. Consequently, the true ground state 
of the system, with free magnetization, at extremely low 
magnetic fields has never been experimentally investi- 
gated. 

The production of Cr BECs in optical dipole traps 
0, 0] allows for the first time to study 5 = 3 spinor 
physics, with a wealth of possible quantum nematic 
phases 111] , depending on the different scattering 
lengths ao.2,4,6 of the molecular potentials S m = 0, 2, 4, 6. 
Contrarily to the cases of Rb and Na where spin depen- 
dent interactions are very small due to similar scatter- 
ing lengths in different molecular channels, Cr has large 
spin dependent contact interactions. One consequence is 
that one can reach the regime where these interactions 
overwhelm the Zeeman effect at experimentally accessi- 
ble magnetic fields (up to 1 mG, compared to typically 
10 /iG for Rb). While large magnetic fields favor a fer- 
romagnetic ground state, below a critical magnetic field 
depending on the interactions the Cr spinor ceases to be 
ferromagnetic [lol which we here observe. 

An important feature of Cr atoms, arising from 
their large electronic spin, is the strength of long- 
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FIG. 1: Chromium BEC spin composition at low fields, revealed 
by Stern-Gcrlach analysis (sec Methods 1). The BEC spontaneously 
depolarizes as the magnetic field is lowered. Absorption pictures af- 
ter 155 ms of hold time, in a field of: a) 1 mG; b) 0.5 mg c) 0.25 
mG and d) mG. At the lowest magnetic field, the spin distribu- 
tion is {17.5 ± 9, 18 ± 4, 14 ± 1.5, 15 ± 3, 17 ± 3, 12.5 ± 4, 6 ± 2} % cor- 
responding to a magnetization of -0.5 gs^B- 



ranged, anisotropic dipole-dipole interactions between 
them. These interactions are too weak to greatly mo dify 
the phase diagram of Cr at low magnetic fields [HI Il2j . 
but they do crucially modify spinor properties. They 
introduce magnetization-changing collisions [l3j |. allow- 
ing the study of spinor physics with free magnetization, 
with an intriguing coupling between the spin degrees 
of freedom and mechanical rotation that resembles the 
Einstein-de-Haas effect [III 14- ltij. They are also ex- 



pected to lead to spin textures, as investigated in [17j . 
similar to domain formation in ferro-magnets. Here, we 
study how dipolar interactions induce magnetization dy- 
namics, as the magnetic field is quenched below an ex- 
tremely low value, which corresponds to a phase transi- 
tion separating two spinor phases of different magnetiza- 
tion. 

In this article, we show that a Cr BEC in which atoms 
are polarized in the lowest single particle energy state 
771 5 = — 3 spontaneously depolarizes below a critical field 
of about 0.4 mG. We observe that the magnetic field be- 
low which depolarization occurs is larger if the peak local 
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density of the BEC is increased by loading it into deep 
2D optical lattices. This is a consequence of an increased 
meanfield energy, which raises the critical magnetic field 
up to 1.2 mG. We also analyze depolarization dynam- 
ics with and without the lattice, and the role of dipolc- 
dipole interactions; we find that depolarization is slower 
in the lattice although the peak density is then higher. 
This counterintuitive behavior is a consequence of a re- 
duction of the average density in the lattice and of the 
non-locality of the mean-field due to long-ranged dipolar 
interactions. In addition, without the lattices, we observe 
an increase of the thermal fraction, as expected from 
non-interacting spinor thermodynamics theory; however 
the observed depolarization of the BEC itself is a feature 
unique to the non trivial phase diagram at low magnetic 
field due to spin-dependent contact interactions. 

We produce Cr BECs by performing evaporative cool- 
ing of atoms in the single particle lowest energy state 
ms = —3, in an optical dipole trap Q. The magnetic 
field during evaporation is sufficiently large that any two- 
body inelastic process is energetically forbidden. After a 
BEC has been obtained in mj = —3, with typically 20 
000 atoms, no discernible thermal fraction, and a peak 
density of 3.5 x 10 14 cm~ 3 in an harmonic trap of fre- 
quencies (320,400,550) Hz, we suddenly (within a 1/e 
time of 8 ms, set by the inductance of the coils and Eddy 
currents) reduce the magnetic field B from 20 mG to an 
extremely low value, typically below 1 mG. 

The magnetic field is controlled by three orthogonal 
pairs of coils, and calibrated by rf spectroscopy. Ex- 
tremely low magnetic fields correspond to magnetic reso- 
nance frequencies of the order of 1 kHz or less. To reduce 
technical noise, and provide better long term field stabil- 
ity, we have implemented an active stabilization of the 
ambient magnetic field (see Methods 1). With active sta- 
bilization, typical shot to shot magnetic field fluctuations 
are 100 yC. 

If the final magnetic field is sufficiently small (below 
0.5 mG typically), we observe a spontaneous depolariza- 
tion of the BEC, while the total number of atoms re- 
mains constant. We measure the population in each of 
the Zeeman sublevels by performing a Stern-Gerlach ex- 
periment (see Methods 2). We finely tune the current 
sent to three orthogonal pairs of coils to maximize depo- 
larization, which in practice precisely pinpoints B = 0. 
Typical depolarization results are displayed on Fig. [T] 

We have repeated the same experiment for a BEC adi- 
abatically (within 15 ms) loaded in the ground state band 
of a 2D optical lattice [18] , superimposed onto the optical 
dipole trap. Our lattice configuration is described in fl9j . 
The lattice depth in each direction is between 25 and 30 
En, where En = 2? ^ A2 is the recoil energy, with A = 532 
nm, and m the atom mass. The vibrational trapping 
frequency in a lattice site is 120 kHz, much larger than 
the chemical potential (jJ,/h =11 kHz); consequently, the 
motion is frozen in two dimensions. We consider that the 
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FIG. 2: Demagnetization curve for 3D BECs and for ID quantum 
gases. The ms — — 3 final population is plotted as a function of the 
magnetic field. Full circles give the results for 3D-BECs and triangles 
for ID tubes. Inset: final gas magnetization as a function of the mag- 
netic field in the lattice, for two orthogonal orientations of the field 
with respect to axis of the lattice tubes. Lines arc gaussian fits to the 
data. 



BEC is split into an array of ID quantum gases, with an 
estimated peak density of 2 x 10 15 cm~ 3 . 

The fraction of atoms remaining in ms = —3 after a 
hold time of 155 ms is displayed on Fig. [2] for the two 
experimental configurations, as a function of magnetic 
field. While we observe strong depolarization at zero 
magnetic field for both configurations, the width of the 
depolarization curve is significantly larger for the BEC 
loaded in the 2D optical lattice than without the lattice. 
A gaussian fit to the experimental data gives a 1/e width 
of 0.4 mG for the BEC, and 1.2 mG for the BEC loaded 
in the lattice. 

We interpret the increased width of the depolarization 
curve in the lattice as an effect of an increase of the mean- 
field interactions. Indeed, the chemical potential of the 
BEC is n/h — 3.8 kHz, and reaches 11 kHz when loaded 
into the lattice. The widths of the depolarization curves 
are therefore roughly proportional to the chemical poten- 
tial. 

At T = 0, and at finite magnetic field B, sponta- 
neous depolarization from ms — —3 becomes energeti- 
cally possible when the energy cost to go from ms = —3 
to ms = — 2 is compensated for by a reduction of in- 
teraction energy between the atoms. In the case of Cr, 
because 04 < a§, this happens below a threshold B c [icij ]: 



gs^BB c « 0.7 



(1) 



with gs ~ 2 the Lande factor of Cr atoms, [ib the Bohr 
magneton and n the density. B c corresponds to a critical 
magnetic field separating the ferromagnetic phase to ei- 
ther a polar phase (spin state (a, 0, 0, 0, 0, 0, /?)), a cyclic 
phase (spin state (a, 0, 0, 0, 0, j3, 0)), or a phase in which 
all Zeeman states are pop ulated, depending on the un- 
known value of a 11]. Phase separation is also pos- 
sible [2(3] • Let us compare the numerical values for B c to 
the values at which we do observe depolarization: for the 
BEC and the lattice respectively, B c calculated from eq. 
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is respectively 0.25 mG and 1.15 mG, which is com- 
parable to the observed 0.4 mG and 1.2 mG 1/e widths. 
The slight discrepancy between experiment and theory 
is consistent with our magnetic field fluctuations. We 
interpret the spontaneous depolarization as due to the 
phase transition at B c . Our work therefore differs from 
|21[ , where depolarization of a thermal Cr gas occurred 
because the thermal energy was converted into Zeeman 
energy through dipolar collisions between atoms. 

We emphasize that technical noise can not mimic the 
increase of the width of the depolarization curve as the 
atoms are loaded in the lattice. Indeed, as the Zeeman 
effect is purely linear for Cr, rf spectroscopy, and more 
generally, the coupling of atoms to time-dependent mag- 
netic fields, are independent of mean-field interactions: 
in a linear spin system, a magnetic field does not couple 
two different molecular potentials. We have also verified 
by rf spectroscopy that the quadratic effect related to the 
tensorial light shift of Cr [22j in the lattice is not respon- 
sible for the increased width of the depolarization curve 

We have investigated the dynamics of depolarization, 
with and without the lattice, as shown on Fig. [3] At 
t=0, the magnetic field control is suddenly switched to 0. 
The magnetic field B(t) does not instantaneously reach 
its final value, and in practice, for the first 50 ms, the 
atoms remain completely polarized as B(t < 50ms) > B c . 
Then, as B(t) < B c , depolarization occurs. Without the 
lattice, depolarization occurs so suddenly (< 5 ms) that 
we can barely resolve it; depolarization is slower in the 
lattice, although the peak density is then much larger. 
As shown below, the typical timescale for depolarization 
is set by the magnetic field resulting from the surround- 
ing dipolar atoms, i.e. by the non-local meanfield due to 



dipolar interactions [14J. For atoms loaded into an op- 



tical lattice, the large increase of the (repulsive) contact 
mean-field forces the cloud to swell in our experimental 
configuration; the overall volume of the cloud is then in- 
creased by a factor of about three, hence reducing the 
dipolar mean-field. A slower depolarization dynamics in 
the lattice is thus a consequence of the non local charac- 
ter of dipole-dipole interactions, and indicates inter-site 
inelastic dipolar couplings in the lattice. 

Another signature for inter-site coupling is provided in 
the inset of Fig. [21 where we show that the depolarization 
is identical for two different orientations 9 of the magnetic 
field, either parallel or perpendicular to the tubes axis. 
In a previous work, where we have studied dipolar relax- 
ation in 2D optical lattices for atoms in ms = 3 at larger 
magnetic fields, the magnetization changing collisions oc- 
cur in a given lattice site, and they are suppressed for 
9 = because of cylindrical symmetry [19j . The situa- 
tion is different here because the experiment is performed 
at much smaller magnetic fields, and dipolar relaxation 
occurs at increasin g in ter-atomic distances for decreas- 
ing magnetic field [13j. At 1 mG, the typical distance 
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FIG. 3: Time evolution of the ms — — 3 population, as a function of 
the time delay between the moment at which the magnetic field control 
is set to zero and the beginning of the Stcrn-Gcrlach procedure, a) 
dynamics for 3D BECs; b) dynamics for ID quantum gases. Full lines 
arc guides to the eye. Depolarization starts after a delay, corresponding 
to the time for B(t) to reach B c (see text). 



at which dipolar relaxation occurs is 300 nm, exceeding 
the lattice periodicity, which breaks cylindrical symme- 
try. A description including inter-sites coupling is hence 
required. 

To account for the depolarization dynamics, we 
consider the time-dependent Gross-Pitaevskii equation, 
starting with a polarized BEC in ms = —3, including 
contact and dipole-dipole interactions. Dipole-dipole in- 
teractions provide the only mechanism to transfer the 
atoms from the maximally polarized ms = —3 state to 
ms = —2. We calculate the dynamics at short times by 
assuming that the population in ms — — 2 remains small, 
and that population in ms > —2 is negligible. Then, 
the coupling between the fields 0j=(— 3,— 2) descri bing the 
ms = (—3,-2) components is given by (see also |2J|): 



,[{x-x')-i{y-y')] {z-z>) 



0-3 (f) 



( 2 ) 

where 7 = -3S 3 / 2 h 2 d 2 /y/2, with d 2 = /^(ffsMs) 2 /4tt, 
and no the magnetic constant. At short times, coupling 
between ms = —3 and ms = —2 atoms is only set by 
dipolar interactions via the term T(r), and |r(r)| 1 is 
the typical time-scale for depolarization at B = 0. This 
timescale is about 3 ms for the BEC case and 10 ms for 
the lattice case, in relatively good agreement with the 
observations (Fig. [3]) at the lowest achievable magnetic 
field. 

We have further monitored the state of the atoms at 
the first stages of depolarization. For this, we have im- 
plemented an alternative Stern-Gerlach procedure (see 
Methods 2), allowing a simultaneous analysis of spin and 
momentum distributions. As shown in Fig. |4j soon after 
depolarization has started the momentum distributions 
of the Zeeman components remain extremely narrow, a 
signature that the system is still Bose condensed. These 
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FIG. 4: Stcrn-Gcrlach monitoring of the spin composition as the 
depolarization develops, for the BEC case. These pictures utilize our 
second Stern-Gcrlach procedure (sec Methods 2) which insures simul- 
taneously spin separation and a measurement of the momentum distri- 
butions of each Zeeman components. The three curves, vertically offset 
for better clarity, correspond to the three absorption pictures shown in 
insert, taken at three steps of depolarization, separated by a few ms. 
Solid lines result from multiple gaussian fits. The width of the nar- 
row peaks shows that, during depolarization, atoms in different spin 
states remain Bosc condensed. A broad additional pedestal reveals the 
presence of a thermal component. 



distributions are narrower than the one of the polarized 
BEC, presumably reflecting the decrease in mean-field 
energy in the depolarized case, since etg is the largest of 
the Cr scattering lengths. We could not detect vortices 
in the ms = —2 component, although they are predicted, 
due to conservation of total angular momentum 11, Tij ]. 
ft is not yet established whether the magnetic field is 
low enough for such vortices to be observed, nor how 
long they should survive in our non cylindrical trap in 
the presence of thermal excitations. This will be further 
investigated in future experiments. 

By measuring the momentum distribution using stan- 
dard time-of-flight techniques without separating the 
Zeeman components, we also observe a decrease of the 
condensate fraction when depolarization is important, al- 
though the temperature (deduced by fitting the thermal 
wings) does not change. Such condensate fraction reduc- 
tion at constant temperature is a direct consequence of 
the release in the constraint of the spin degrees of free- 
dom, which lowers the critical degeneracy temperature 
T c |25[. Although the reduction of T c might be under- 
stood in the framework of non-interacting spinor thermo- 
dynamics, it is worthwhile noting that depolarization of 
the BEC can not be understood in such framework, as 
a non-interacting BEC with free ma gne tization remains 
ferromagnetic for all magnetic fields [26j. 

In conclusion, this work represents the first investiga- 
tion of quantum magnetism where magnetization is free 
and spin-dependent contact interactions set the many- 
body mean-field ground state. We have observed the 
magnetization dynamics of a BEC due to dipolar interac- 
tions, as the magnetic field is quenched through a critical 
value corresponding to (at T=0) a quantum phase transi- 
tion. We have established that demagnetization is set by 



the dipolar mean-field. Demagnetization leads to a non- 
polarized BEC, although the spin state which is reached 
does not necessarily reveal the nature of the ground state, 
due to diabaticity at the phase transition, and to ther- 
mal excitations: temperature (~ 100 nK) is comparable 
to the difference of energy between the spinor phases at 
B = 0. 

We will next try to observe the rotation necessarily 
induced by demagnetization, in the spirit of the Einstein- 



de-Haas effect [11|, |14H16| . We will also investigate how 



this rotating system thermalizes at low magnetic fields, 
to possibly reach quantum nematic (mean-field) phases 
(10| . When the BEC is loaded into the 2D lattice, we 
produce an array of mesoscopic samples (40 atoms per 
site) at low magnetic field, high density, and with large 
spin-dependent interactions; such conditions are known 
to be favorable to the observation of non-classical spin 
states 27] or fragmented BECs 28]. 
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METHODS 



Magnetic field control 

The three components of the magnetic field along or- 
thogonal directions of space are stabilized using an active 
cancellation of the AC and DC magnetic field fluctua- 
tions: three large, 3-turn rectangular coils, with a typ- 
ical size of 1.5 x 1.7 m, and with their center located 
about 1 m away from the BEC, along the 3 spatial direc- 
tions, are used as compensation coils. The magnetic field 
close to the experimental chamber is measured using a 3 
axis fluxgate magnetic sensor (Bartington Mag 03MC). 
The sensor is located outside the experimental chamber, 
about 15 cm away from the BEC. The three components 
of the magnetic field are compared to a computer con- 
trolled set-value, with a resolution of 50 fj,G. The error 
signal is sent to a feedback loop which is a simple propor- 
tional controller, and the current in each coil is regulated 
with a 1 Amp push-pull type voltage-controlled current 
amplifier. 

Without active compensation, the characteristics of 
the magnetic noise are the following: the main noise con- 
tribution comes from the 50 Hz AC noise of magnitudes, 
up to 4 mG peak/peak in one of the directions. This noise 
originates from the different equipment located around 



the experiment. The DC magnetic field fluctuates up 
to 2 mG during the day time. This is much larger than 
the natural fluctuations of the earth magnetic field which 
are of the order of 100 to 200 fj,G. These fluctuations are 
mainly caused by activities in the laboratory (positions 
of metallic objects) and outside the laboratory (elevators, 
cars). 

When the active stabilization is switched on, the AC 
noise is decreased below 100 /iG peak/peak and the DC 
magnetic field value is stabilized to better than 20 /jG at 
the position of the sensor. We have also measured the 
residual magnetic field fluctuations with a second inde- 
pendent sensor located 20 cm away from the first one. 
We observe that the drifts of the DC magnetic field stay 
below 100 /jG over 1 hour. At this position, though, 
the typical AC fluctuations are 500 fiG peak/peak. This 
comes from the fact that some of the AC noise sources 
are not very far from the BEC (around 1 meter) so that 
the AC noise is not perfectly spatially homogeneous on 
a 20 cm scale. We have nevertheless checked by rf spec- 
troscopy that this AC noise is efficiently screened by the 
metallic experimental chamber, so that the residual fluc- 
tuations at the BEC position are given by the 100 /j>G 
DC fluctuations. The experimental spectrum shown in 
Fig. [2] corroborates this estimate. 

Stern and Gerlach procedures 

We have implemented two different Stern-Gerlach pro- 
cedures to monitor the spin populations at low magnetic 
fields. For both procedures, a small (10 mG) magnetic 
field B Q is first adiabatically applied to the atoms, with a 
rise time of 10 ms. We insure that B a is large enough that 
a magnetic field gradient can also then be applied to the 
atoms to separate the Zeeman components adiabatically, 
i.e. without changing the spin composition. 

In the first Stern-Gerlach procedure, the applied mag- 
netic field gradient is small (0.25 G/cm), and the mag- 
netic field felt by each Zeeman component remains al- 
ways small. The Zeeman components then separate very 
slowly. To prevent them from falling and keep them in 
the imaging field of view, their horizontal motion is chan- 
neled in one direction by the horizontal laser beam used 
for producing and trapping the BEC, and only the ver- 
tical beam is removed during the procedure. An absorp- 
tion image is taken after 45 ms of ID expansion in the 
magnetic field gradient, right after the horizontal opti- 
cal trap is switched off. Due to the very small magnetic 
field employed, absorption imaging offers a good absolute 
determination of the number of atoms in each Zeeman 
component, as the Zeeman shifts are small compared to 
the transition linewidth for absorption imaging. We ver- 
ify this by the fact that the sum of the populations in all 
Zeeman components is constant while the gas depolar- 
izes. Unfortunately, the momentum distribution of each 
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Zeeman component is lost when using this procedure, be- 
cause the atoms have then moved out of the focus of the 
imaging system. This first procedure was used for all 
results in this paper, except the ones shown in Fig. 0J 

A second Stern-Gerlach procedure was used for Fig. 2) 
We apply a stronger magnetic field gradient (1 G/cm). 
The Zeeman components separate sufficiently rapidly 
that the optical trap can be released almost simultane- 



ously. The Zeeman components are separated in 5 ms, 
and this faster procedure also gives access to the momen- 
tum distribution of each Zeeman component, although 
the actual determination of their relative atom number 
is not as good as the one in the first procedure, because 
of the presence of a relatively strong magnetic field (a 
few G) when absorption imaging is performed. 



